Ab initio Studies on the Interplay between Spin-Orbit Interaction and Coulomb 

Correlation in Sr2lr04 and Ba2lr04 
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Ab initio analyses of A2lr04 (A=Sr, Ba) are presented. Effective Hubbard-type models for Ir 5d 
t2g manifolds downfolded from the global band structure are solved based on the dynamical mean- 
field theory. The results for A=Sr and Ba correctly reproduce paramagnetic metals undergoing 
continuous transitions to insulators below the Neel temperature Tjv. These compounds are classified 
not into Mott insulators but into Slater insulators. However, the insulating gap opens by a synergy 
of the Neel order and significant band renormalization, which is also manifested by a 2D bad metallic 
behavior in the paramagnetic phase near the quantum criticality. 



PACS numbers: 71.30.+h,71.20.-b,71.15.-m 

-Introduction. Electron correlation effects in quasi- 
two-dimensional systems, especially transition metal 
compounds, have been one of the central issues in con- 
densed matter physics since the discovery of the cuprate 
superconductors. While the effective electron-electron in- 
teraction is large in the 3d transition metal compounds, 
it becomes weaker in compounds built from heavier tran- 
sition metal elements, e.g. Sd, because of their more spa- 
tially extended orbitals. On the other hand, the strength 
of the spin-orbit interaction (SOI) increases and thus in 
5d compounds it is competitive with other characteristic 
energy scales: The interplay between SOI, on-site repul- 
sion, crystal-field splitting and inter-site hopping opens 
a new field of research in strongly correlated compounds. 

Ir compounds offer such model systems. In particular, 
Sr2lr04 and Ba2lr04 are isostructural to La2Cu04 im- 
plying a two-dimensional anisotropy[ll|2|. Since ten-fold 
degenerate Ir bd bands arc partially occupied nominally 
by an odd number (five) electrons, the insulating behav- 
ior of Sr2lr04[3[, with antiferromagnetic (AF) order be- 
low 250Ki^i5|, indicates an important role of electronic 
correlations. Based on the band structure of Sr2lr04, 
Kim et aL[6| have shown that the t2g manifold near the 
Fermi level is well separated from the eg one. SOI lifts 
the t2g degeneracy splitting it into four, mostly filled, de- 
generate states indexed with the angular momentum j ~ 
3/2, and a doubly degenerate j — 1/2 state filled with 
nearly one electron. They also suggested that low-energy 
excitations are accounted for by a single-band Hubbard 
model with two pseudospin states jz = ±1/2 expressed 
by a linear combination of the three t2g states; \j = 
1/2; jz = ±1/2) = {\yz, ±a) Ti\zx, ±ct) t \xv, T<j))/^, 
where a describes the spin. They proposed that Sr2lr04 
provides a realization of (weakly canted) AF Mott in- 
sulator similar to the parent materials of the high- Tc 



cuprates. Recent resonant X-ray scattering measurement 
also supports this scenario [7[. 

Ba2lr04 is expected to have a narrower bandwidth be- 
cause of the larger ionic radius of Ba. However, the AF 
transition around 240K is similar to the Sr compoundj^. 

In this Letter, we present a detailed ah initio analysis of 
A2lr04 (A=Sr, Ba). To predict accurately the electronic 
structure from the first principles, we employ a recently 
proposed three-stage scheme j8l-llCl|. First, we obtain the 
global band structure using the density functional the- 
ory (D FT) within the generalized gradient approximation 
(GGA)[ll|. Second, using a Wannier projection on the 
Ir t2g target bands we eliminate the states far from the 
Fermi level. We employ the constrained random phase 
approximation (cRPA) [9| to obtain the screened interac- 
tion parameters for the downfolded six-orbital Hubbard- 
type model. Third, we solve the derived low-energy ab 
initio model and study its dynamics with the dynami- 
cal mean- field theory (DMFT)[l3]. This general ah ini- 
tio framework may be called Multi-energy-scale Ah initio 
scheme for Correlated Electrons (MACE). MACE repre- 
sents a general and hierarchical framework [SJ, which en- 
ables ah initio calculations of strongly correlated electron 
systems. Among possible choices of MACE, here, we em- 
ploy the combination of DFT and DMFT. 

Our results reveal that, in agreement with experi- 
ments, Sr2lr04 and Ba2lr04 show a similar behavior, 
both with a continuous transition from a low temperature 
(T) AF insulator to a high-T paramagnetic meiaZ[l3j, 
but with a smaller energy scale (bandwidth, interaction 
strength) for Ba2lr04. Despite not being concomitant, 
the gap opening is induced by the AF transition, in con- 
trast to La2Cu04 but possibly similar to Nd2Cu04[lJ]. 
Therefore, we classify the iridates as Slater insulators. 
The effective single-band and square-lattice behavior 



makes these materials valuable reference systems for the 
cuprate superconductors. Indeed, we find their paramag- 
netic phases to be strongly renormalized 2D metals. Our 
results also provide insight into the proposed possible 
high Tc superconductivity [l5| and giant magnetoelectric 
effects [la| in these compounds. 

-ab-initio downfolding. In Fig. [l] we show the band 
structure of Sr2lr04 and Ba2lr04 obtained by the den- 
sity functional calculation. We used the PBE exchange- 
correlation functional [uj and the augmented plane wave 
and local orbital (APW-l-lo) method including the 
spin-orbit coupling as implemented in the WIEN2K 
program [13, El. The muffin tin radii (Rmt) of 2.23, 1.97, 
1.74 bohr for Sr, Ir and O were used for Sr2lr04 and 2.34, 
1.92, and 1.70 bohr for Ba, Ir and O for Ba2lr04, respec- 
tively. The maximum modulus for the reciprocal vectors 
^max was chosen such that i?MT-f''max = 7.0 and a 10 x 
10 X 10 fc-mesh in the first Brillouin zone was used. 

The crystal structure of Ref.[19[ was used for Sr2lr04. 
The obtained band structure is similar to previous 
studiesja, |2Q|. The experimental lattice parameters, a 
and c, of Ba2lr04 arc 3% and 4% larger than those of 
Sr2lr04[2], while the Ir-O-Ir angle is yet to be deter- 
mined experimentally. Therefore, we have optimized the 
structure requiring any atomic force to be less than 0.2 
mRy/bohr. The obtained Ir-O-Ir angle is 164°. The dif- 
ference between the optimized (153°) and experimental 
(159°) Ir-O-Ir angles for Sr2lr04 suggests a few degrees 
overestimate of the Ir-O-Ir in Ba2lr04 as well. 

We then constructed the Wannier functions for the 
i2g-like bands, using the WIEN2Wannier 
wannier90 J22l | codes. In Fig. [l] we supermipose 



and the 
the 

Wannier-interpolated bands on the original GGA bands. 
The hopping between the d^y orbitals is well represented 
by a two dimensional tight-binding model while the dyz 
and dzx bands are quite one-dimensional. The spin-orbit 
coupling mixes the cubic harmonics to form eigenstatcs 
of the pseudospin j. We used the j, jz basis, which diag- 
onalizes the onsite part of the Hamiltonian and in which 
the local Green function is to a good approximation di- 
agonal, in the subsequent calculations. The hopping am- 
plitudes are listed in the Supplemental material. In the 
right panels of Fig. [1] we show the partial density of 
states (pDOS) of the j = 1/2 and j = 3/2 states. The 
center of gravity of the j = 1/2 bands is higher than 
that of j = 3/2 bands, and its contribution is dominant 
around the Fermi level. It is interesting to note that 
two one-dimensional bands (dyz and dzx) and one two- 
dimensional band (dxy) are mixed up equally and make 
quite two-dimensional j = 1/2 bands. The width of the 
j = 1/2 band is about 1.5 eV, a few times smaller than 
that of the dx^-.y2 band of the cuprates. If we replace 
Sr with Ba, the unit cell expands, and the band width 
becomes smaller. 

Next, we evaluate the interaction parameters by 
cRPaQ. We used the Density Response Code (DRC)^ 
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FIG. 1: (Color online) (Left) GGA band structure (solid 
lines) together with the Wannier interpolated band of the t2g 
states ((red) dotted lines) for Sr2lr04 (for the experimen- 
tal structure (a) and optimized structure (b)) and Ba2lr04 
(c). (Right) Partial density of states for j = 1/2 ((red) solid 
line), j = 3/2, j^ = ±l/2((blue) thin line) and j = 3/2, j^ = 
±3/2((black) dotted line) states. 



recently developed for the Elk branch of the original 
Exciting FP-LAPW code[25|. Below we sketch the main 
steps of the calculation. DRC utilizes the time dependent 
DFT to compute density susceptibility, x(r, r',a;) = 
X^'S'(r,r',a;) + //dridrs x^'^lr, ri, w)(w(ri, rz) -I- 
/'''=(ri,r2,a;))x(r2,r',w). Here w(ri,r2) = |ri - r2|"^ is 
the bare Coulomb interaction and /^"^ is a dynamical 
exchange-correlation kernel ignored in the conventional 
cRPA. 



,KSt 



X"-'\r,r',uj) = Y^ 
j-j' 



'-^{fj-Ir)^ri^)rr{v')U^')r,{v) 



UJ — [eji — Ej) + ^0^ 



(1) 
is the susceptibility of the noninteracting Kohn-Sham 



TABLE I: Effective on-site Coulomb ([/) /exchange (J) inter- 
actions between two electrons on the same Ir site in the t2g 
model for all the combinations of Ir-5d orbitals (in eV). 
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electrons. Here, e^ and fj are the energy and oc- 
cupancy of the eigen-state V'ji ^-nd X) runs over all 
pairs of bands but excludes the cases of j^j' both be- 
longing to the target t2g bands. Then the screened 
Coulomb interaction defined by M^(r,r',w) = u(r,r') + 
JJ dridr2 w(r,ri)x(ri, r2, w)w(r2, r', w) yields Hubbard- 
U parameters ("4- index [/" matrix): 

Ua^^s = ^im dridr2 w*^{ri)w*p{r2)W{ri, r2, uj) X 

X w^{ri)'ws{r2) (2) 

with Greek letters representing a combined index of 
band and translation. We used the band-disentanglement 
methodj2(| to unhybridize the t2g bands from the rest. 

Here, to save the numerical cost, we constructed the 
simplified crystal structure with one formula unit per unit 
cell moving the inplane oxygen atoms to the symmetric 
positions. We took 100 unoccupied bands and 5 x 5 x 5 fc- 
and <7-meshes in the calculation of the dielectric function. 
The double Fourier transform of x(r, r', w) was done with 
the |G -I- <7| = 3.5 [1/a.u.] cutoff which gives ^540 and 
~470 G- vectors for Ba2lr04 and Sr2lr04, respectively. 

In Table U we list the Coulomb and exchange interac- 
tions {Uij = Uijij and Jij = Uijji, where i and j denote 
the t2g orbital index) in the cubic harmonics basis. From 
those we have deduced an approximate parametrization 
in terms of Slater parameter Fq, F2, while fixing the 
ratio F4/F2 to the atomic value of 0.625. We have 
found Fa = 1.933(1.623) eV and F2 == 2.266(2.068) 
eV well parametrize the interaction matrix for Sr2lr04 
(Ba2lr04). From these Slater parameters we have cal- 
culated the pairwise interactions Uij in the j, j^ basis 
(see the Supplemental material). The expanded unit cell 
for Ba2lr04 makes the crystal field splitting between the 
t2g and Bg smaller and the screening becomes more effi- 
cient. This nearly cancels the band narrowing, leaving, 
the ratio of the band width to the interaction strength 
similar upon replacing Sr with Ba. This similarity leads 
to similar behaviors of the two compounds contrary to 
naive expectation. 

-DMFT analysis. In DMFT, the inter-atomic corre- 
lations are neglected, and the system is mapped self- 
consistently on an effective impurity model, which we 



solve with a Monte Carlo method[27| keeping only the 
density-density terms of the on-site interaction. We have 
performed two types of calculations. First, we have in- 
vestigated the Mott insulator scenario. For this purpose 
the system was constrained to the paramagnetic phase 
and the interaction was treated as a free parameter of 
simplified form Uij = U (this approximation will be jus- 
tified a posteriori). Second, we studied the real mate- 
rials with the interaction strengths at their respective 
first-principles cRPA values (including the orbital de- 
pendence). We allowed symmetry breaking towards a 
long-range magnetic order and studied its evolution as a 
function of temperature. For all the interaction strengths 
studied, we have found that the both materials behave 
effectively as a one- band systems with the j = 1/2 be- 
ing half-filled, while the j = 3/2 bands completely filled. 
Therefore the j = 1/2 intra-band interaction is the only 
parameter, which controls the charge gap opening. 



(a)Sr IrO (experimental) 



(b)Ba IrO (optimized) 
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FIG. 2: (Color online) Non-magnetic LDA+DMFT result for 
Sr2lr04 (a) and Ba2lr04 (b) at T = 1/40 eV. 



One-particle spectra in the paramagnetic phase ob- 
tained for various interaction strengths are shown in 
Fig. [51 The metal-insulator transition occurs ai U ^ 
2.3(1.8) eV for Sr2lr04 (Ba2lr04) while the intra-j=l/2 
repulsion from cRPA amounts to 1.96 (1.65) eV for 
Sr2lr04 (Ba2lr04). Thus for the cRPA interaction, both 
compounds remain metallic if paramagnetic. These met- 
als, however, have large mass renormalization factors 
3.5(6) for Sr2lr04 (Ba2lr04) at U = [/crpa and 
1/80 eV. Here, a rough lower bound of Z is given by 
1 — ImS(iwi)/a;i at the lowest Matsubara frequency 



Z 
T 
Z 

UJl. 

Next, we fix the interaction to the cRPA values. We 
drop the paramagnetic constraint and allow the system 
to stabilize an AF order, without altering the Bravais 
lattice to keep the calculation feasible. Such an in-plane 



checkerboard order is consistent with the experiment, but 
the inter layer stacking configuration differs [7|. In Fig. |3l 
we show the result for the one-particle spectra. At high 
temperatures we obtained metallic paramagnetic solu- 
tions, while both of Sr and Ba compounds undergo a 
transition to the AF phase at around 810 (690) K for 
the Sr2lr04 (Ba2lr04). As the staggered magnetiza- 
tion grows below T/v, a gap continuously opens in the 
one-particle spectra. The continuous transition with a 
nearby quantum critical point and the resultant enhanced 
AF fluctuations near it naturally cause the renormalized 
large- .Z metals described above. This is in contrast with 
the first-order nature claimed in the recent variational 
Monte-Carlo calculation [28|. While the calculations pre- 



(a)Sr IrO (experimental) 



(b)Ba2lrO^ (optimized) 
-■ — I — <- 



diet the two compounds to order at similar temperatures, 
the actual values are overestimated by a factor 2-3. This 
is not surprising given the lack of the inter-site fluctua- 
tions in DMFT and the layered nature of the materials. 

The gap opening can be understood considering the 
quasiparticle equation e* — eq ~ ^'ai^a) = 0; where e* is 
the quasiparticle energy, eo is the bare energy, cr = ± is 
the pseudospin index and E' is the real part of the self- 
energy. S' consists of a frequency independent Hartree 
term and frequency dependent 'wiggle' at low frequency, 
responsible for the mass renormalization in PM phase. 
In the traditional Slater picture only the Hartree part of 
S' is considered and the gap opens when the 'E'_^_ — S'_ 
splitting exceeds the bare bandwidth. Here, the Hartree 
splitting, about 30% smaller than the bare bandwidth, 
is not enough to open the gap. The necessary ingre- 
dient is the low- frequency 'wiggle' in S^(cj) and its a- 
dependent shift towards positive (negative) frequencies 
for the empty (fllled) orbitals. Similar cr-dependent shift 
in the case of an AF Mott insulator was observed in [29[ . 
While the Hartree splitting is a simple consequence of a 
static polarization and does not depend on particular or- 
dering, this cr-dcpendcnt shift of S'g.(a;) is characteristic 
of the AF order, where the empty and filled orbitals on 
the neighboring sites have the same a, generating a hy- 
bridization repulsion between the empty and filled bands. 
The large mass renormalization Z > 3.5 — 6 present al- 
ready in the PM phase cooperatively helps the gap open- 
ing. In this regard, we note that Z is larger than Z ~ 2 
for SrVOa and CaVOs, typical 3d correlated metals [30[. 
While DMFT is not expected to predict Tjv accurately, 
it correctly captures the formation of the one-particle 
gap induced by the developed AF correlations, which in 
real systems may exist even in some temperature interval 
above T^. On the other hand, in the absence of SOI, the 
system is metallic in the LDA+U calculationlg. Given 
the fact that LY)A+U usually overestimates the insulat- 
ing gap, we should indeed have a metallic solution in 
LDA+DMFT, if we neglect SOI. 

-Conclusion. The electronic structures of Sr2lr04 
and Ba2lr04 calculated from first principles combined 
with the DMFT show that they have AF and insulat- 
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FIG. 3: (Color online) AF LDA-f DMFT result for Sr2lr04 (a) 
and Ba2lr04 (b). The Hubbard U estimated by cRPA (1.96 
for (a) and 1.6 for (b), respectively) is used. The overestimate 
of 0.05 < Tn < 0.1 in eV in comparison to the experimental 
values ~ 250K may be due to the neglect of the spatial fluc- 
tuations in DMFT. Staggered magnetization as a function of 
temperature is plotted in (c). 



ing ground states. These iridates undergo continuous 
transitions to paramagnetic metals above T/v, indicat- 
ing that they are to be classified as Slater insulators, 
in agreement with the available experiments. Though 
they have an essential difference from the Mott insula- 
tor, they are not simple Slater insulators either, because 
strongly renormalized bad metals emerge in the para- 
magnetic phase. The present Slater insulators are the 
consequence of substantial cooperation of Mott-type cor- 
relation effects. The strongly-renormalized paramagnetic 
metal adjacent to the AF Slater insulator opens a pos- 
sibility of unexplored correlation effects under the inter- 
play of the spin-orbit interaction. The similarity and dis- 
similarity to the cuprates elucidated here offers intrigu- 
ing reference systems to understand the superconducting 
mechanism in the cuprates when carriers arc doped. 
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